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Previously  a  lab  scale  catalytic  autothermal  reformer  (ATR)  capable  of  operating  at  pressures  from  6  to 
50  bar  was  constructed  and  tested.  The  objective  of  the  experimental  program  was  to  maximize  H2  pro¬ 
duction  per  mole  of  02  supplied  (H2(OUt)/02(in)).  In  this  companion  paper  a  1-D,  heterogeneous,  numerical 
model  is  developed  and  tested  for  simulating  the  high  pressure  ATR.  The  effects  of  molar  steam  to  carbon 
(S/C)  and  oxygen  to  carbon  (0 2/C)  ratios  are  studied  and  optimal  operating  conditions  are  identified  for 
three  system  operating  pressures;  6,  28  and  50  bar.  Experimental  optimal  conditions  and  model  results 
are  compared  and  found  to  be  in  close  agreement.  The  optimal  conditions,  however,  predicted  by  the 
model  at  pressures  of  28  and  50  bar  have  higher  S/C  ratios  and  produce  higher  H2(OUt)/02(in)  yields  than 
the  experimentally  determined  optimums.  A  sensitivity  analysis  consisting  of  9  model  parameters  is  also 
performed.  The  model  is  most  sensitive  to  the  activation  energy  of  the  two  steam  reforming  reactions 
used  in  the  model  and  the  operating  parameter  02/C. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Significant  research  and  development  effort  has  been  expended 
to  advance  hydrogen  as  the  energy  carrier  of  tomorrow.  While  it 
is  likely  that  future  hydrogen  will  be  produced  from  renewable 
resources  such  as  biomass  or  electrolysis  powered  by  geother¬ 
mal,  photovoltaic  or  wind,  currently  the  most  economic  way  to 
produce  hydrogen  is  by  catalytic  thermochemical  reforming  of 
hydrocarbons.  Methane  in  particular  has  been  the  primary  source 
for  hydrogen  production  and  techniques  developed  for  methane 
reforming  can  aid  in  transitioning  to  sustainable  renewable  sources 
of  hydrogen  for  tomorrow.  Over  the  last  few  decades,  kinetic 
modeling  of  methane  reforming  through  various  thermochemical 
techniques  has  been  reported  in  the  literature.  Results  have  identi¬ 
fied  methods  to  increase  efficiency  and  provided  insight  useful  in 
improving  reactor  design. 

In  light  of  advancing  reforming  techniques,  there  has  been 
recent  interest  in  constructing  small  scale,  underwater,  fuel  reform¬ 
ing  stations.  Applications  include  powering  underwater  sensory 
equipment  and  recharging  unmanned  submersibles.  In  concept, 
these  remote  reforming  stations,  using  onboard  fuel  cells  to  gen¬ 
erate  power,  would  be  deployed  on  the  ocean  floor  where  near 
natural  gas  seeps  or  methane  hydrates.  Methane  from  these  sources 
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would  be  collected  and  purified  to  remove  catalyst-poisoning  H2S 
and  then  reformed  to  produce  a  hydrogen  rich  gas  for  the  fuel  cell. 
While  methane  would  be  available  in  abundance,  a  major  limitation 
for  operating  a  reformer  and  fuel  cell  on  the  seafloor  is  having  oxy¬ 
gen  available  at  sufficiently  high  concentrations.  Additionally,  since 
the  unit  would  be  deployed  on  the  ocean  floor  and  operate  at  ambi¬ 
ent  pressure,  design  of  components  would  need  to  accommodate 
the  effects  of  elevated  pressure  on  the  physical  system  and  reformer 
and  fuel  cell  chemistry.  An  alternative  to  operating  the  reformer  at 
ambient  pressure  would  be  to  operate  the  reformer  in  a  submerged 
pressure  hull  that  would  allow  the  system  to  operate  at  a  lower  pres¬ 
sure.  The  latter  system  would  require  that  methane  at  the  ambient 
pressure  of  the  ocean  floor  be  decompressed  and  that  unwanted 
byproducts  be  recompressed  for  removal  from  the  pressure  hull. 
Operating  the  reformer  at  ambient  pressure  was  selected  over  this 
alternative.  For  this  fuel-rich,  oxygen  limited,  reforming  applica¬ 
tion  it  is  of  particular  interest  to  maximize  moles  of  H2  produced 
per  mole  of  02  supplied. 

Results  of  an  experimental,  laboratory-investigation  of  methane 
reforming  at  simulated  sea  floor  pressure  conditions  were  reported 
in  a  previous  paper  by  the  same  authors  [1].  This  paper  reports  a 
companion  modeling  effort  that  explores  elevated  pressure  reform¬ 
ing  via  numerical  simulation. 

I  A.  Modeling 

This  section  briefly  describes  modeling  work  which  has  been 
conducted  on  the  three  main  thermochemical  catalytic  reform¬ 
ing  techniques;  steam  reforming  (SR),  partial  oxidation  reforming 
(POX)  and  autothermal  reforming  (ATR).  It  should  be  noted  at  the 


554 


M.A.  Reese  et  al.  /  Journal  of  Power  Sources  195  (2010)  553-558 


Nomenclature 

Q  concentration  of  gas  species  i  (mol  m-3 ) 

C?n  concentration  of  inlet  gas  species  i  (mol  m-3 ) 

c^g  average  specific  heat  of  gas  (J  kg-1  K-1 ) 
dp  catalyst  pellet  diameter  (m) 

Ddpi  dispersion  coefficient  of  gas  species  i 

Di  gas  diffusivity  of  species  i  (m2  s-1 ) 
h  heat  transfer  coefficient  between  catalyst  bed  and 

gas  phase  (W  nrr3  K-1 ) 
i  species  index 

j  reaction  index 

A Hj  heat  of  reaction  j  (J  mol-1 ) 

mave  average  molar  mass  of  gas  into  which  gas  species  i 
is  diffusing  (gmol-1) 
mi  mass  of  gas  species  i  (g  mol-1 ) 

P  system  pressure  (bar) 

Pi  partial  pressure  of  gas  species  i  (bar) 

r2-  reaction  rate  of  gas  species  i  (mol  kg-1  catalyst  s-1 ) 

Rj  reaction  rate  of  reaction  j  (mol  kg-1  catalyst  s-1 ) 

Sh  heat  transfer  area  per  unit  volume  of  catalyst  bed 

(m-2  m-3) 
t  time  (s) 

T  temperature  of  gas  phase  (I<) 

Vn  gas  inlet  temperature  (K) 

Ts  temperature  of  solid  phase  (K) 

u  superficial  gas  velocity  (m  s-1 ) 

z  coordinate  (m) 

Greek  symbols 

£  void  fraction  of  the  catalyst  bed 

A  effective  thermal  conductivity  of  catalyst  bed 

(Wnrr1 1<-1) 

Peat  density  of  catalyst  bed  (kg  catalyst  m-3 ) 

~pl  average  gas  density  (kg  m-3 ) 

rbed  tortuosity  of  the  bed 

Subscripts 

bed  catalyst  bed 

cat  catalyst 

8  gas 

i  gas  species  (H2,  02,  CO,  C02,  H20,  CH4) 

j  reaction  index  ( 1  -4) 


outset  that  SR  and  POX  have  been  implemented  in  industry  for 
some  time  [2].  As  such,  much  literature  can  be  found  describing 
the  thermodynamics  and  kinetics  of  relevant  reactions.  In  compar¬ 
ison,  relatively  few  papers  exist  on  the  kinetics  of  ATR.  The  following 
is  a  brief  sample  of  thermochemical  reforming  models  which  have 
been  developed  and  studied. 


1.2.  Steam  reforming 

In  1989,  Xu  and  Froment  developed  a  rate  equation  and  kinetic 
data  for  methane  steam  reforming  (MSR)  [3].  Their  study  was  con¬ 
ducted  in  the  temperature  range  of  573-823  K  with  gauge  pressure 
varying  from  0.3  to  lObar.  A  conventional  tube  reactor  was  used 
and  Ni  based  catalyst  employed.  Their  work  formed  a  basis  for 
much  modeling  work  in  the  MSR  field  and  has  also  carried  over 
into  POX  and  ATR  work.  Hou  and  Hughes  [4]  developed  slightly 
different  rate  equations  and  kinetic  data  and  tested  their  model  at 
1.2-3.0  bar.  Both  SR  models  consisted  of  three  global  reactions:  two 


endothermic  steam  reforming  reactions: 

CH4  +  H2O  -*  CO  +  3H2,  AH1(298)=  206 kj mol"1  (1) 

CH4  +  2H20  ->•  C02+4H2,  AH2(298)=  165 kj mor1  (2) 

and  the  exothermic  water  gas  shift  reaction: 

CO  +  H20  C02+H2,  AH3{298)=  -41  kjmor1  (3) 

1.3.  Partial  oxidation  reforming 

Modeling  of  POX  was  outlined  by  De  Groote  and  Froment  in  1996 
[5].  Their  model  was  based  on  the  use  of  a  Ni  catalyst  and  consisted 
of  rate  expressions  for  nine  reactions  including  Eqs.  (l)-(3)  and 
an  additional  six  reactions  describing  methane  combustion,  C02 
reforming,  carbon  deposition,  methane  cracking,  and  carbon  oxi¬ 
dation.  The  reaction  set  was  employed  in  a  one-dimensional,  fixed 
bed,  adiabatic,  heterogeneous  model  based  on  work  from  Froment 
and  Bishoff  [6]. 

Other  models,  such  as  the  POX  membrane  reactor  proposed  by 
Basile  et  al.  in  2001  [7],  used  only  reactions  in  Eqs.  (l)-(3)  and  the 
combustion  of  methane: 

CH4  +  202^  C02+2H20,  AH4(298)  =  —802  kj  mol-1  (4) 

Like  many,  the  model  proposed  by  Basile  et  al.  was  one-dimensional 
and  valid  near  atmospheric  pressure  (1.2  bar) 

1.4.  Autothermal  reforming 

Chan  and  Wang  [8]  provided  an  overview  of  ATR  performance 
using  chemical  equilibrium  models  and  identified  operating  con¬ 
ditions  to  maximize  H2  yield  while  minimizing  CO  production  and 
avoiding  carbon  deposition.  To  further  the  research  in  this  field, 
Hoang  and  Chan  described  ATR  kinetics  in  2004  [9].  This  would  be 
the  first  of  their  series  of  papers  on  kinetic  modeling  of  MSR,  POX 
and  ATR  [10-12].  Their  reaction  scheme  for  ATR  was  a  combination 
of  SR,  based  on  Xu  and  Froment  [3],  and  POX  reforming  reactions. 
Including  only  those  reactions  with  significant  rates  [12]  led  to  a 
reduced  reaction  set  based  on  Eqs.  (l)-(4).  All  reforming  papers 
by  Hoang  et  al.  used  a  two-dimensional,  heterogeneous  model  at 
atmospheric  pressure. 

In  2006,  Lee  et  al.  [  13  ]  presented  results  from  a  one-dimensional, 
pseudo-homogenous  model  based  on  Hoang  and  Chan’s  reaction 
scheme.  The  effects  of  steam  to  carbon  (S/C)  and  oxygen  to  carbon 
(02/C)  molar  ratios  in  the  reactant  streams  to  the  reformer  were 
studied  at  atmospheric  pressure. 

De  Groote  and  Froment  reported  simulating  a  25  bar  ATR  con¬ 
dition  while  validating  their  POX  model,  but  this  was  only  a  single 
run  and  no  study  involving  pressure  effects  was  presented  [5]. 

In  1998,  Olsvik  and  Hansen  presented  a  paper  on  high  pres¬ 
sure  autothermal  reforming.  Their  work,  however,  dealt  primarily 
with  fluid  flow,  combustion  kinetics,  and  soot  formation.  No  data 
on  product  gas  composition  for  the  various  reformer  conditions  of 
interest  were  presented  [14]. 

Therefore,  while  many  papers  dealing  with  thermodynamics 
and  kinetics  of  various  reforming  techniques  exist,  very  few  papers 
have  been  found  which  deal  with  high  pressure  methane  ATR. 

1.5.  Objectives 

The  objectives  of  this  study  were  to  adopt  a  suitable  model  for 
high  pressure  ATR  in  order  to  compare  modeling  results  with  exper¬ 
imental  results  previously  reported  [1  ]  and  to  conduct  a  sensitivity 
analysis  of  the  model. 
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Fig.  1.  Schematic  layout  of  1  -D  reactor. 

2.  Materials  and  methods 

A  one-dimensional  ATR  model  based  on  mass,  energy,  and 
species  balances  was  developed  which  included  kinetic  rate  expres¬ 
sions  for  chemical  reactions.  The  model  was  solved  using  the 
software  package,  COMSOL  Multiphysics  3.4  (COMSOL,  Inc.,  Los 
Angeles,  CA)  based  on  the  finite  element  method  (FEM). 

2.1.  Model  description 

The  model  was  based  on  reactions  in  Eqs.  (l)-(4)  and  included 
seven  gas  species:  methane,  oxygen,  steam,  carbon  monoxide,  car¬ 
bon  dioxide,  hydrogen,  and  nitrogen.  A  schematic  layout  of  the 
1-D  reactor  model  is  shown  in  Fig.  1.  Tz  represents  the  gas  phase 
temperature,  Ts,z  represents  the  solid  phase  temperature,  and  Qz 
represents  the  concentration  of  gas  species  i  at  a  given  location  z.  L 
represents  the  length  of  the  catalyst  bed  (0.61  m)  and  the  nodes  at 
which  the  solution  is  calculated  are  represented  by  dots. 

2.2.  Model  assumptions 

Model  development  included  the  following  assumptions:  (1) 
convective  transport  dominates  diffusive  transport  in  the  axial 
direction  for  thermal  and  mass  diffusion  (Peclet  number  >140  and 
>40,  respectively)  allowing  the  diffusion  terms  to  be  omitted,  (2) 
the  thermal  conductivity  of  the  gas  phase  is  much  smaller  than  that 
of  the  solid  phase  catalyst  bed  and  can  be  omitted,  (3)  radial  varia¬ 
tions  are  considered  negligible  and  therefore  all  radial  components 
are  omitted. 

2.3.  Model  equations 

The  kinetic  rate  equations  and  kinetic  parameters,  equilibrium 
constants  and  adsorption  constants  were  obtained  from  a  publica¬ 
tion  by  Eloang  et  al.  that  employed  Engelhard’s  Ni-0309S  catalyst 
[11].  With  the  assumptions  mentioned  above  the  three  governing 
equations  that  describe  the  model  are  listed  below.  Eqs.  (5)-(7) 
are  the  gas  phase  species  balance,  the  energy  balance  for  the  solid 
phase,  and  the  energy  balance  for  the  gas  phase,  respectively. 

=  ~u~^  +  Pcatt/,  i  =  1  to  number  of  species  (5) 

4 

PbCpcat^  =  A.0  +PcatJ2(-AHj)Rj+Shh(T-Ts)  (6) 

j=i 

_ dT  _ dT  _  . 

£Pgcpg~fo  —  ~uPgcpg~^  +Shh(Ts  -  T )  (7) 

where 

Ddpi  =  e(^+ 0.5  d„uj  (8) 

and 

Tbed  —  (9) 

D;  =  1  X  10s  ■  r1-75  .  (( Vmgve)  +  (Vmi))1/2  ^  (10) 

P((20.1)1/3  +  (P;  *  20.1)1/3) 


2.4.  Boundary  conditions 

The  following  boundary  conditions  were  specified  in  the  model: 
at  the  reactor  inlet,  z  =  0  T  =  Tin;Ts  =  Tm;Q  =  Cjn  (11 ) 

dT  dT  dC~ 

and  at  the  reactor  outlet,  z  =  L  —  =  0;  A.— ^  =  0;  -1  =  0  (12) 

dz  dz  dz 

2.5.  Initial  conditions 

The  model  was  initialized  with  the  following  values.  The  inlet 
temperature  of  the  gas  species,  Tn,  was  set  to  the  experimental 
inlet  temperature  of  270  °C  [1].  The  catalyst  bed  temperature  was 
given  a  linear  profile  from  270  °C  at  the  reactor  inlet  to  700  °C  at 
the  reactor  outlet.  The  CH4  flow  rate  (2.51  min-1),  N2  flow  rate 
(0.5 1  min-1 ),  appropriate  steam  to  carbon  (S/C)  and  oxygen  to  car¬ 
bon  (O 2/C)  (variable)  ratios  and  pressure  (6, 28,  or  50  bar)  were  also 
specified.  While  the  inlet  temperature  and  gas  flow  rates  were  held 
at  their  initial  values  over  the  course  of  a  run,  the  catalyst  bed  tem¬ 
perature  profile  was  not  and  the  profile  changed  significantly  as  the 
reactions  occurred. 

2.6.  Solver  and  mesh  size 

A  time  dependent  solver  was  used  to  propagate  the  model 
solution  from  the  initial  conditions.  Using  the  temperature  pro¬ 
file  assigned  to  the  catalyst  bed  the  reaction  started  at  the  outlet 
of  the  reactor  and  propagated  toward  the  inlet,  until  the  reaction 
zone  stabilized.  The  results  from  this  apparent  steady  state  solution 
determined  with  the  time  dependent  solver  were  saved  and  used 
as  the  initial  conditions  for  the  stationary  solver.  This  was  done  to 
verify  the  time  dependent  solution.  If  the  results  converged  using 
the  stationary  solver,  a  steady  state  solution  for  the  set  of  operating 
condition  was  assumed  to  have  been  reached.  A  parametric  solver 
was  subsequently  used  to  determine  solutions  for  different  values 
of  model  parameters. 

The  mesh  size  was  critical  to  ensure  accurate  results  from  the 
model.  If  the  mesh  was  too  large,  steep  gradients  were  missed  and 
faulty  calculations  made.  Mesh  refinement  was  performed  to  deter¬ 
mine  if  the  mesh  size  in  use  was  small  enough. 

3.  Results  and  discussion 

3.1.  Comparison  of  modeling  and  experimental  results 

The  results  of  the  ATR  model  were  compared  to  the  experimental 
results  previously  acquired  [  1  ].  The  experimental  optimal  condition 
for  maximizing  H2(OUt)/02(in)  at  6  bar  (02/C  =  0.44;  S/C  =  3.35),  28  bar 
(02/C  =  0.48;  S/C  =  3.0)  and  50  bar  (02/C  =  0.48 ;  S/C  =  3.0)  for  exper¬ 
imental  and  model  output  gas  composition  is  compared  on  the 
parity  plot  shown  in  Fig.  2.  All  data  points  are  moles  of  gas  species 
out  divided  by  moles  of  gas  species  in  unless  otherwise  noted.  Under 
these  conditions,  experimentally,  one  mole  of  02  produced  4.16, 
3.43  and  2.84  moles  of  EI2  at  6, 28  and  50  bar,  respectively.  The  mod¬ 
els’  correlation  to  experimental  results  for  these  conditions  is  high. 
In  fact  the  model  was  able  to  predict  the  decreased  performance 
at  increased  pressure  very  well:  all  predicted  results  were  within 
±10%  of  experimental  results  except  the  moles  of  CO(OUt)/CH4(in). 
It  should  be  noted  that  Boudouard’s  equation  was  not  included 
in  the  model  and  therefore  the  model  is  not  suitable  for  predict¬ 
ing  carbon  formation.  This  approach  was  taken  since  the  operating 
conditions  tested  in  previous  experimental  work  [1]  did  not  favor 
carbon  formation. 

Since  the  model  predicted  outlet  gas  composition  in  high  corre¬ 
lation  to  the  experimental  results,  the  model  was  used  to  look  at  a 
broader  range  of  operating  conditions  than  experimentally  tested. 
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Fig.  2.  Parity  plot  of  experimental  results  and  kinetic  model  results  at  6,  28  and 
50  bar.  (All  data  points  are  moles  of  gas  species  out  divided  by  moles  of  gas  species 
in  unless  otherwise  noted.) 


Figs.  3-5  show  the  moles  H2(OUt)/02(in)  over  the  full  range  of  con¬ 
ditions  modeled  for  6,  28,  and  50  bar,  respectively.  Only  moles  of 
H2(0ut)/02(in)  is  shown  since  this  was  the  focus  of  the  accompanying 
experimental  work. 

The  experimental  work  identified  reformer  operating  points 
that  were  not  stable  and  resulted  in  the  high  temperature  reaction 
zone  being  blown  out  of  the  reactor.  This  behavior  was  observed  in 
the  model  results  as  well.  For  a  fixed  pressure  and  0 2/C  ratio  the 


Fig.  4.  Moles  of  H2(out)/02(in)  model  results  for  28  bar  (various  0 2/C  and  S/C  ratios). 


model  would  attempt  to  parametrically  solve  a  series  of  different 
S/C  ratio’s.  With  a  starting  value  of  2.0  the  S/C  ratio  would  be  gradu¬ 
ally  increased  with  a  small  step  size  to  a  final  value  of  4.5.  In  nearly 
every  test,  however,  at  a  particular  S/C  ratio,  the  high  temperature 
reaction  zone  would  move  downstream  through  the  reactor  with 
large  steps  when  the  S/C  ratio  was  increased  ever  so  slightly.  When 
this  occurred,  modeling  for  the  given  02/C  ratio  was  stopped  and 
then  restarted  using  the  next  02/C  ratio  to  be  tested  with  a  S/C 
value  of  2.0.  For  example  in  Fig.  3  at  an  0 2/C  ratio  of  0.35  when  the 
S/C  ratio  reached  2.6  this  phenomenon  was  observed.  At  that  point 
modeling  was  restarted  at  an  O 2/C  ratio  of  0.40  and  a  S/C  ratio  of 
2.0. 

It  was  observed  that  as  pressure  increased  a  higher  S/C  ratio  was 
attainable  for  a  given  O 2/C  ratio.  This  was  due  to  both  the  higher 
resulting  reactor  temperature  and  the  lower  resulting  space  veloc¬ 
ity  which  helped  keep  the  flame  front  from  blowing  out  the  back 
end  of  the  reactor. 

The  optimal  yield  and  operating  condition  for  moles  of 
H2(out)/°2(in)  predicted  by  the  model  was  4.09  at  6  bar  (02/C  =  0.45; 
S/C  =  3.1),  3.65  at  28  bar  (02/C  =  0.5;  S/C  =  3.7)  and  3.08  at  50 bar 
(02/C  =  0.50;  S/C  =  4.4).  At  6  bar  the  predicted  optimum  and  the 
corresponding  02 /C  and  S/C  ratios  are  very  comparable  to  the  exper¬ 
imental  values.  At  28  and  50  bar  the  predicted  optimum  was  higher 
and  also  occurred  at  a  higher  S/C  but  comparable  O 2/C  ratio  when 
compared  to  experimental  values.  This  suggests  that  the  model  was 
able  to  run  under  conditions  having  greater  S/C  ratios  than  exper¬ 
imentally  observed  at  elevated  pressures.  Being  able  to  run  with 
a  higher  S/C  ratio  meant  higher  molar  FI2(OUt)/02(in)  yield’s  were 
attainable. 

Fig.  6  shows  the  predicted  moles  of  H2/02(in)  along  the  length  of 
the  reactor  at  various  S/C  ratios  an  O 2/C  ratio  of  0.5  and  a  pressure 
of  50  bar.  The  model  predicted  a  very  steep  increase  in  moles  of 
H2/02(in)  near  the  inlet  with  a  much  more  gradual  increase  towards 
the  backend  of  the  reactor.  The  slope  of  the  moles  of  H2/02(in) 
curve  at  the  front  of  the  reactor  is  determined  by  two  factors:  the 
cooling  effect  of  steam  which  is  injected  at  270  °C  and  the  space 
velocity  which  increased  with  increasing  S/C  ratios.  Therefore,  with 
each  increase  in  the  S/C  ratio  the  moles  of  H2/02(in)  curve  became 
more  gradual  at  the  front  of  the  reactor  in  spite  of  reaching  higher 
final  values.  The  most  extreme  example  is  observed  at  a  S/C  ratio 
of  4.5  where  the  curve  is  clearly  pushed  downstream.  Also  note, 
although  pressure  effects  are  not  shown  in  Fig.  6,  the  difference  in 
the  FI2/02(in)  yield  between  successive  S/C  ratios  decreased  as  pres¬ 
sure  increased.  This  suggests  the  model  is  more  sensitive  to  the  S/C 
ratio  at  low  pressures  than  it  is  at  elevated  pressures. 

Figs.  7-9  shows  the  dry  percent  volume  of  each  species  in  the 
reformate  gas  over  a  range  of  S/C  ratios  at  an  O 2/C  ratio  of  0.5  and 


M.A.  Reese  et  al.  /  Journal  of  Power  Sources  195(2010)  553-558 


557 


Fig.  6.  Moles  of  H2/02(in)  model  results  for  50  bar  along  the  reactor  length  (various 
S/C  ratios). 


pressures  of  6,  28  and  50  bar.  The  trends  at  each  pressure  are  the 
same:  an  increased  S/C  ratio  led  to  increased  percentages  of  H2  and 
CO2  and  corresponding  decreases  in  CH4,  CO,  and  N2.  As  pressure 
increased,  however,  the  percentages  of  CO  and  C02  remained  nearly 
identical  while  the  concentrations  of  CH4  and  N2  increased  and  H2 
decreased. 

3.2.  Sensitivity  analysis 

Figs.  10  and  11  illustrate  the  sensitivity  of  the  kinetic  model  with 
respect  to  9  model  parameters:  the  steam  to  carbon  ratio  (S/C),  oxy¬ 


Fig.8.  Kinetic  model  results:  dry  gas  concentration  (various  S/C;  02/C  =  0.5;  28  bar). 


Fig.  9.  Kinetic  model  results:  dry  gas  concentration  (various  S/C;  02/C  =  0.5;  50  bar). 


Fig.  10.  Moles  of  H2(0Ut)  sensitivity  analysis  (28  bar). 


gen  to  carbon  ratio  (0  2/C),  reaction  rate  of  reaction  1  (RR1 ),  reaction 
rate  of  reaction  2  (RR2),  activation  energy  of  reaction  1  (El ),  activa¬ 
tion  energy  of  reaction  2  (E2),  overall  heat  loss  (HL),  conductivity 
of  the  catalyst  bed  (Cat  k),  and  pressure  (P).  The  reaction  rate  and 
activation  energy  of  reaction  3  and  4  had  a  negligible  effect  in  all 
studies  and  therefore  were  not  presented.  Five  different  model  out¬ 
puts  (moles  of  H2,  CO,  C02,  CH4,  and  %  H2)  were  studied  by  varying 
the  9  model  parameters  ±25%  from  their  standard  value  when  run 
at  28  bar  a  S/C  =  3.0  and  an  0 2/C  =  0.48. 

The  moles  of  H2(0Ut)  was  most  sensitive  to  E2  followed  by  02/C 
and  El.  A  25%  increase  in  E2  and  0 2/C  resulted  in  a  60%  decrease 
and  >15%  increase  in  moles  of  H2(0Ut),  respectively.  Similarly  a  25% 
decrease  in  RR1  decreased  the  moles  of  H2(0Ut)  by  over  15%.  The 
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sensitivity  analysis  of  the  %  H2(0Ut)  as  well  as  the  moles  of  CO(OUt) 
and  C02(0ut)  followed  the  same  trend  as  H2(0Ut).  The  trends  in 
the  sensitivity  analysis  for  the  moles  of  CH4(0Ut)  were  nearly  the 
opposite  of  the  H2(0Ut)  trends  and  can  be  seen  in  Fig.  11. 

The  moles  of  CH4(0Ut)  was  most  sensitive  to  E2  and  0 2/C  followed 
by  El.  By  increasing  E2  25%  the  moles  of  CH4(0Ut)  increased  by  over 
80%,  by  increasing  0 2/C  25%  the  moles  of  CH4(0Ut)  decreased  by  over 
40%  and  finally  by  decreasing  El  25%  the  moles  of  CH4(out)  decreased 
by  nearly  30%. 

In  summary  the  model  output  parameters  were  most  sensitive 
to  E2,  02/C,  and  El,  respectively.  In  all  cases,  increasing  E2  had  the 
opposite  effect  of  increasing  El  and  02/C.  Therefore  increasing  E2, 
or  decreasing  El  and  02/C,  resulted  in  an  increase  in  the  moles  of 
CH4(0Ut)  and  a  decrease  in  the  moles  of  H2(0Ut),  CO(OUt),  and  C02(OUt) 
and  %  H2(0Ut).  The  model  tended  to  be  the  least  sensitive  to  HL  and 
RR2. 

4.  Conclusions 

A 1-D,  heterogeneous,  kinetic  model  was  adapted  for  high  pres¬ 
sure  ATR  and  tested.  Experimental  optimal  conditions  and  model 
results  at  6,  28  and  50  bar  were  compared  and  found  to  be  in 
close  agreement.  At  elevated  pressures  the  model  predicted  stable 
reformer  operation  at  S/C  ratios  higher  than  were  experimentally 
possible  for  a  given  02/C  ratio.  Optimal  conditions  predicted  by 
the  model  at  pressures  of  28  and  50  bar  had  higher  S/C  ratios  and 
produced  higher  H2(OUt)/02(in)  yields  than  the  experimentally  deter¬ 
mined  optimums.  The  optimum  conditions  predicted  by  the  model 
at  6  bar  was  4.09  moles  of  H2(OUt)/02(in)  (02/C  =  0.45;  S/C  =  3.1),  at 
28  bar  was  3.65  moles  of  H2(out)/02(in)  (O2/C  =  0.5;  S/C  =  3.7)  and 
at  50  bar  was  3.08  moles  of  H2(-out)/02(in)  (02/C  =  0.50;  S/C  =  4.4). 


A  sensitivity  analysis  performed  on  the  model  showed  that  the 
model  was  most  sensitive  to  the  activation  energy  of  the  two  steam 
reforming  reactions  used  in  the  model  and  the  operating  param¬ 
eter  02/C.  By  using  modeling  techniques  the  effects  of  individual 
system  variables  on  system  performance  were  studied  and  opti¬ 
mal  operating  conditions  were  identified  for  three  system  operating 
pressures. 
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